Ordering of anisotropic polarizable polymer chains on the full many-body level 
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We study the effect of dielectric anisotropy of polymers on their equilibrium ordering within 
mean-field theory but with a formalism that takes into account the full n-body nature of van der 
Waals forces. Dielectric anisotropy within polymers is to be expected as the electronic properties 
of the polymer will typically be different along the polymer than across its cross section. It is 
therefore physically intuitive that larger charge fluctuations can be induced along the chain than 
perpendicular to it. We show that this dielectric anisotropy leads to n-body interactions which 
can induce an isotropic-nematic transition. The two body and three body components of the full 
van der Waals interaction are extracted and it is shown how the two body term behaves like the 
phenomenological self-aligning-pairwise nematic interaction. At the three body interaction level we 
see that the nematic phase that is energetically favorable is discotic, however on the full n-body 
interaction level we find that the normal axial nematic phase is always the stable ordered phase. 
The n-body nature of our approach also shows that the key parameter driving the nematic-isotropic 
transition is the bare persistence length of the polymer chain. 



I. INTRODUCTION 



Semiflexible polymers, just like rod-like anisotropic nematogens, form orientationally ordered liquid crystalline 
mesophases, spanning the concentration regime of chiral nematics f5] and all the way to line hexatics Among the 
most important examples of semi flexible polymer ordering is the double stranded DNA molecule [U [5] which makes 
a plethora of liquid crystalline phases in aqueous solutions. These ordered phases of DNA aqueous dispersions are 
pertinent to the construction of DNA nanoassemblies [B] and have recently found a promising new application in the 
form of virus-like particles |7j. 

A rigorous theory of nematic ordering of long rigid nematogens in dilute solutions was first proposed in a seminal 
work by Onsager |8j. His approach was later adapted to semiflexible polymers by Semenov and Khokhlov [9j and 
within a different formal background by Ronca and Yoon |TT]. A fleld-theoretical reformulation of the nematic- 
isotropic transition for semi flexible polymers was then pursued by Gupta and Edwards [T as well as by Tkachenko 
and Rabin |12) . These approaches are all based on the assumption of pairwise additivity of the interaction potential, 
assumed to be of a general nematic form stemming from the Onsager steric interaction. For short range interactions 
this is of course a reasonable assumption. However, it is well known that electrostatic as well as van der Waals (vdW) 
interactions usually engender long range non-pairwise additive interactions that in general do not conform to this 
approximation |13) . 

Furthermore, polymers in general show an anisotropic dielectric response, so that the dielectric functions along the 
axis and perpendicular to the axis differ. This is nicely seen e.g. in the case of DNA [Hj, where amplitude of the 
dielectric function in the axial and the radial directions are very different because of the strong optical anisotropy of 
DNA |15j . This dielectric anisotropy has important consequences as vdW interactions between anisotropic media lead 
in general to torques induced by electromagnetic field fluctuations as first realized by Weiss and Parsegian [T71 [TB] . 
The details of this effect are complicated but it persists not only in the non-retarded but also in the retarded limit 

[IMH. 

Existence of anisotropic vdW interactions now poses a question about their possible role in ordering transitions. 
While the fundamental importance of the vdW force for the gas-liquid transition is well recognized we now investigate 
the possible role of anisotropic vdW interactions for orienational ordering transitions. As it is not clear off-hand 
whether the effect is signiflcant, it is of utmost importance that the many-body pairwise non- additive nature of 
vdW interactions be fully recognized in the theoretical framework. We thus formulate the theory of orientational 
ordering of semi flexible polymers due to anisotropic vdW interactions in the way that allows for a complete and 
exact resummation of vdW interaction energy to all levels, including pairs, triplets, quadruplets etc. Whatever the 
conclusions of this theory, they can certainly not be undermined by the improper treatment of the vdW effect. 

In order to assess the effect of anisotropic vdW interactions we analyze a rather simple case of semiflexible polymer 
ordering driven exclusively by these interactions. We thus study the effects beyond the usual Onsager steric anisotropy 
ansatz. However, the most fundamental difference with the previous work is that we calculate exactly the contribution 



2 



of vdW interactions to all orders within the mean-field approximation. Usually vdW interactions are treated on a 
pairwise level, or possibly three body Axilrod- Teller level [22 , but never more. We propose an analytical approach 
that allows us to sum the van der Waals interactions to all orders and solve the corresponding mean-field theory 
exactly. It should be noted that the long range nature of van der Waals interactions implies that our mean-field 
treatment should be reasonably accurate, especially in dense systems. Pursuing this approach we show indeed that 
the anisotropic vdW interactions can engender an orientational ordering transition in a polymer solution of stiff chains. 



II. ELECTROSTATIC FIELD ENERGY 



We start with a polymer described as a continuous chain within the arclength parametrization as x(s) so that the 
tangent vector is defined as 

t(s)=i(s). 

The chain is assumed to be composed of an anisotropic macroscopic dielectric material that has a dielectric response 
function £|| and ej^ parallel and perpendicular to the chain, respectively. The electrostatic energy of the polymer- 
solvent system is thus described as 



£es = he 



(V0(x))^dx+ie|| / i,ij V,(/)(x)Vj-<?!)(x)rfx+ ie_L / (<5.j - V,(/.(x)Vj-0(x)dx. (I) 



Assuming now that the polymer is much thinner, of radius a, then its length we get the final expression of the 
electrostatic energy as 



£es = / (V0(x))' dx + ^ / ds{{e^-e^) V0(x(s))2 + (e|| - e^) t,t, V,0(x(s)) V,0(x(5))) . (2) 

Here the first integral is now over the whole space and the second one is a line integral along the polymer chain. We 
introduce the scalar polymer density field pijx.) and the tensor polymer nematic field (Tij(x) defined as 

p(x)= / ds5(x-x(s)), (3) 

and 

a,,{^)= I dsS{^-xis))U{x{s))t,{^{s)), (4) 

where L is the total length of the polymer. We can now write the electrostatic energy as 

£bs[0(x)] = i / e~y(x)V,0(x)Vj0(x), (5) 
Jv 

where the composite dielectric function is now given by 

Sij{x) = + 7ra^(e_L - £w) p(x)) dij + 7ra^(£|| - e_L) aij{x). (6) 
In thermodynamic equilibrium electric field fluctuations give rise to the partition function 

= J 2?[0(x)] exp(-^f£5['/'(x)]). (7) 

Since the integral over the electrostatic field, in the limit where the polymer fields are both independent and constant 
i.e. p(x) = p and cr^ (x) = aij , is of a general Gaussian type it can be evaluated explicitly giving for the corresponding 
free energy 

f d^k 
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This is the form of the free energy corresponding to thermal fluctuations of the electrostatic field that we will use later. 
It is equivalent to the zero Matsubara frequency or static van der Waals interactions as is well known j22| . In the full 
Lifshitz theory [IB] with assumed complete magnetic homogeneity of the space, the above formula is generalized to 



n=0 



where the n summation (the prime indicates that n = term has a weight of 1/2) is over the imaginary Matsubara 
frequencies ^„ = 2TTnkBT/h, where is the Boltzman constant, T is the absolute temperature, and h is the Planck 
constant divided by 27r. 

III. POLYMER ENERGY 

The polymer chain itself is not featureless and will be modeled within the framework of the semiflexible Kratky- 
Porod model where the polymer conformational energy is given by 



£kp ~ \K I ds t(s)^ with the constraint t(s)'^ = 1, 



(10) 



where the last constraint stems from the fixed length of the monomers. Instead of dealing with this constraint locally, 
which is a difficult task, we simplify the model as suggested by Edwards and Gupta [Ij, by taking it into account 
globally in the form 

{t{sf) = 1. (11) 
This then leads to the effective partition function of the polymer chain in the form 

Skp^\K i{sfds + \\ {t{sf-l)ds, (12) 



where A is a Lagrange multiplier to be evaluated self-consistenly from the global constraint of Eq. ( 11 ). The polymer 
partition function is then given by the form 

Zkp^ j V[^{s)] aqy{-p£KpHs)]) with j V[^{s)] ±{sf eiqi{-P£Kp[^{s)]) = Zkp- (13) 



IV. CONSTRAINT "ENERGY" 



In order to evaluate the total free energy of the system we need to insert the proper constraints into the total energy 
of the fluctuating electrostatic field and polymer configurations. That the two are coupled is evident from Eq. (|6|. 
In the final partition function one thus needs to insert the following delta functions that can be represented with the 
functional delta transforms 



j V[p{Ti)] 5 ^/9(x) - jds (5(x - x(s)) j = j 2?[p(x)]X>[r(x)] exp J dx r(x)p(x 



)-i Ids r(x(s)) y (14) 



where r(x) is the corresponding auxiliary field, and 



Jv[a{x)] S |^f7,,(x) - Jds <5(x - x(s))t,(x(s))tj(x(5)) 



I?[(t(x)]I?[s(x)] exp ^ Jdx s,y(x)CTij(x) - i Jds Sij{x{s))U{x{s))tj{x{s))j , (15) 

where the auxiliary field is now the tensor s(x) with components Sjj(x) and we adopt the Einstein summation 
convention. When using these constraints in the compact definition of the complete partition function we will also 
perform a Wick rotation in the complex plane for the two auxiliary fields as well as redefine both fields in such a 
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way that they scale with the inverse temperature. Since they are arbitrary fields anyway this does not introduce any 
inconsistencies. 

These two constraints can then be rewritten in a form of an effective " constraint partition function" of the form 

Zco= J P[p(x)]P[r(x)]2?[a(x)]P[s(x)]exp(-/3fco[r(x),p(x),s,,(x),ai,(x)])), (16) 



where 



£co['"(x),p(x),s(x),c7(x)] = - / dxr(x)p(x)-i / rfx Sjj(x)a-ij(x) + 

Jv Jv 

+ Ids r(x(5)) + i Ids Sij{^{s))ti{^{s))tj{^{s)). (17) 
Jo Jo 



V. TOTAL PARTITION FUNCTION 

In view of the above the total partition function can be cast into the following form 



Zt = J I?[(/.(x)]I?[x(s)] J I?[p(x)]I?[r(x)]I?[a(x)]I?[s(x)]exp(-/3fT), (18) 



where 



i / £.,(x)V.0(x)V,</.(x) + 
Jv 

pL pL pL pL 

+ \K / y.{sfds+\X / {±{sf-l)ds+\ ds Sij{x{s))ti{x{s))tj{yi{s)) + ds r{x{s)) 
Jo Jo Jo Jo 

- / dx r(x)p(x) - i / dx Sij(x)c7jj(x). (19) 
Jv Jv 

We now assume that both the polymer density field and the polymer orientational field are constant in space, which 
affectively means that we introduce a mean-field ansatz. In that case we already evaluated the functional integral 
over ^{x). We next tackle the polymer part being a functional integral over the polymer configurational field x(s). 
In order to do that we first of all introduce the Rouse or polymer Fourier decomposition as 

x(a) = / ^'^(^)e-'"^ (20) 
The complete polymer part of the partition function can thus be written as 

Z^(s) = exp (i/3AL - ^rL) JJ / dx{cj) exp (-/3fpo(x(a;))) (21) 

where L is the total length of the polymer chain, i.e. L = Nb and 

£po{x{u))) = ^aij{ui) Xi{uj)x* j{uj), (22) 

with 

aij{uj) = u? {{Ku? + \)5ij + Sij) . (23) 

The functional integral over the Rouse modes is now Gaussian and can be evaluated explicitly yielding up to irrelevant 
multiplicative constants 



{det aij{Lo)) 

The intermediate partition function has thus been derived in the form 



^x(«) = . .a/, • (24) 



Zt = J V[p{x)]V[r{x)]V[a{x)]V[s{x)] exp (-^Jv) , (25) 
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where the intermediate free energy entering the above definition can be written as 

PPt = J log (i.jhkj) + \L j logdeta,j (a;) -PV rp- PV \s,ja,j + prL - \l3\L. 



(26) 



In order to derive this form we had to assume that all the relevant fields, i.e. r,p,Sij and atj, are constant in space. 
Taking into account the connection between the monomer density p and the length of the polymer chain, i.e. 

L Nb 

^^^^nb^p, (27) 

we derive the intermediate free energy density as 

l3J^T{ra>,s,(T) ^ ^ j -^^\og{eij{p,aij)kikj) + \nb j '^\ogdeiaij{uj, s^^)- (3 {rp + \sij(Jij)+l3 {r - \\) nb. (28) 
with 

£jj(p, cr) = (e™ + 7ra^(ej_ - e^) p) (5ij + 7ra^(e|| - ej_) (Tij and aij{u:, s) = uj^ {{Kuj'^ + \)5ij + Sij) . (29) 
As indicated, this still has to be evaluated in the thermodynamic equilibrium as a function of all four polymer fields. 

VI. MEAN-FIELD EQUATIONS 

The second integral can be evaluated easily in an explicit form in terms of the eigenvalues of s and of ct, 
yielding for the free energy density 



The minimization is now carried out with respect to all the fields, i.e. s, cr, p, r, X. Notice that the first term in the 
above is the term due to the interactions between the polymer chain with itself, however the following terms are 
just the free energy for a polymer which is suitably constrained. The minimization with respect to the constraint 
parameters s^, r and A gives 



These equations can now be used to eliminate the variables and p, the resulting free energy being independent of 



r. The first equation in Eq. (31 1 along with that of Eq. (32) implies that a must obey the constraint 



E] o"^ = Tr (7 = p = nb. (33) 

The resulting effective free energy is thus a function solely of a and some simple algebra gives 

P^rjcr) ^ pPvdw{(y) {nbYkBT _^ 

V V 32 K ^ ^ ' 

The zero Matsubara frequency van der Waals part of the free energy can be written up to terms independent of a 

as 

^'^^^ = ^ y ^27r"^3 [1'^ i^w + 7ra^(ej. - £«;) nb) + ln(fc^ + aaijhkj)] (35) 
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with 



a. 



7ra^(e|| — ej_) 



(36) 



and we have used the relation p = nb. 

It now remains to compute the zero Matsubara frequency van der Waals free energy which means we must compute 
an integral of the form 



W{B) = J d^k \n{Bijkikj). 



To proceed we notice that 



where 



dW{B) 

dBpq 



f 

Jo 



A3 



/„„ = / dk- 



■^ij ki kj 



(37) 



(38) 



(39) 



Here k denotes the unit vector and A is a ultraviolet cut-off corresponding to length scales below which the electro- 
magnetic field fluctuations are cut off. We now note that 

rf3kexp(--k2)-^ = / k''dkeM-^k')Ip<iiBij) = J-I,,{B,j). (40) 

and use the identity 

/I dt f 1 

d^keM-^^^)kpkg -eM-^Bijkikj) = -{2w)i dt [I + tB)-gUet{l +tB)-i , (41) 

where / denotes the identity matrix. Putting all of this together we obtain 



Ipg{B) = 27r / dt{l + tB)-^ det(/ + tB)-'^ 
Jo 



dW{B) 

dBj,a 



Thus we derive the effective potential W{B) in the form of 



W{B = J + ao-) = 47r / 
Jo 



dt 
~t 



-det(/ + tB)-5 



= 47r 



dt 
~t 



- det{I + t{I + aa))~^ 



(42) 



, (43) 



Using this result the zero Matsubara frequency van der Waals part of the free energy can then be written in a 
compact form as 



P^vdw{(T) A3 



V 127r2 
so that the total free energy is obtained as 



PC 

In (e^ -I- 7ra^(e_L - e^)n6) -)- / 

Jo 



dt 



1 



t \{l + t)i 



-det[I + t{I + aa)]~^ 



(44) 



V 



127r2 



In (sy, + 'na?{e^ - £w)nb) + — ( 3- - det[/ -|- t{I + aa)]~^ ] 

Jo t \(l+t)2 ) 



The zero Matsubara frequency van dor Waals free energy may now be expanded in terms of a, the expansion up to 
order o? corresponding to the pairwise approximation for the van der Waals interaction. Carrying out this expansion 
to third order in a yields up to terms independent of u, 



/3.FT(a) A3 



V 



{nbfkBT 
32 K 



"Tr CT - ^(2 Tr + (Tr of) + ;^(8 Tr ^3 + 12 Tr Tr a + (Tr of) + 0{a^) 
3 30 315 



+ 



+ 



Tr 



(46) 
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It is interesting to compare our free energy functional at the order of o? with that derived within the same mean-field 
approximation for nematic interactions between polymer where the pairwise interaction energy is given by [I], 



5" / / dsds'{t{s) X t(s'))^'5(x(s) - x(s')). (47) 







The form of this interaction makes it energetically favorable for the polymers to align locally in a parallel fashion. In 
the mean field approximation this gives a cr dependent mean-field free energy of nematic interaction 

y = 2^"^*- ~ Tr cr ] + ^ Tr a , (48) 

the part coming from the constrained polymer partition function being identical as a function of a in the two cases. 
Recalling that as Tr cr = nb is fixed we can identify an effe ctiv e nematic interaction parameter for the van der Waals 



interaction by identifying the coefficients of Tr a in Eqs. (461 and (48) to obtain 



Ue = — ^ — ■ (49) 

907r2 ^ ' 

The effective nematic interaction is thus proportional to the square of the relative anisotropy polarizability a and 
thus crucially depends on the non-isotropic nature of the polymer material. Note also that the sign of the interaction 
is always positive. This ties up with the fact that polarizable cylinders have a preference to align in a parallel fashion, 
thus favoring nematic order. However only stiff molecules with non-isotropic dielectric response, such as DNA or 
carbon nanotubes, can be expected to exhibit this particular tendency toward nematic ordering. 

The variational free energy density per unit volume /(cr) = J^t{<^)/V, can be rewritten making a change of variables 
a ^ nb a' where now the constraint is Tr cr' = 1 and the effective free energy per unit volume is obtained up to terms 
independent of cr as 



°° dt f 1 

T 



- dct[/ + t{I + 7ct')]"5 + c Tr o- 



(50) 



The dimensionless parameters 7 and c are defined as 



, ^(g|| , S n^nb) 

J = nba = and c=-— — — -, (51) 

Sro + iy{e\\ - e^) 8 A^lbb 

Here v — mra^b is the volume fraction of the polymer and Ip = K /bhsT is the dimensionless bare persistence length 
of the polymer in units of the monomer length b. From the definition it is clear that < 7 < (ey — e_L)/£||. This simple 
fact is a crucial point arising from our analysis. The fact that the effective interaction parameter saturates at high 
volume concentrations means that a phase transition cannot be induced by simply increasing the concentration of 
polymer as in the case of models with self-aligning-pairwise interactions T. Physically this result can be interpreted 
by the fact that n-body interactions tend to frustrate the ordering which is induced by the pairwise term. 
In terms of the volume fraction we can write the variable c as 

'""^ (52) 



SA^a'^blp 



Furthermore we assume that the monomer size b and the radial size of the polymer a are of the same size. The cut-off 
A in connection with vdW interactions is associated with the breakdown of the continuum dielectric description and 
is usually taken to be of the atomic size [32]. We assume the scaling form A ~ 27rA/a, where A is the ratio between 
the atomic and the macro-molecular sizes A and a. This means that we can write 

^ " C. (53) 



647r2A3 Up 647r2A3 Ip 



and we notice that the precise values of the parameters a and b can be absorbed into the redefinition of the persistence 
length of the polymer, i.e. Ip — > lp(b/a). In the analysis that follows the dimensionless parameter c is expected to 
be small as by definition v < 1 and Zf, > 1. 
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VII. NEMATIC TRANSITION 



The parameter 7 is small in the limit where e|| w ej_ and/or when the volume fraction ly is small. However for strong 
contrasts in ey and e_L it can be of order 1. It is convenient to measure the components of the polymer's dielectric 
tensor in terms of the dielectric constants of the solvent, we write en = e^eu and eu = £^£11, which now gives 



7 



eii - 1 



1 



1 



l + i/(e|| -1) 



(54) 



Therefore 7 is a monotonic function of and the absolute value of 7 increases upon increasing v. In order to proceed 
we write the tensor a as 



where Q is the order parameter tensor 



Q 




(55) 



(56) 



Here S is the standard uniaxial nematic scalar order parameter and T is the measure of the biaxality. 

At the pairwise level, the van der Waals interaction to order can be mapped onto a nematic interaction Eq. 



(49), a uniaxial nematic alignment of the polymers is favoured and indeed Gupta and Edwards [Ij found in their 
analysis that only these phases are thermodynamically stable. Thus in what follows we will examine the possibility 
of uniaxial ordering to all orders in a and thus set T = in the ansatz Eq. ( 56 ) . Note that in this case we must 
have —1/2 < 5* < 1. While within this ansatz positive S corresponds to a normal nematic phase, the lower branch 
with negative S implies excess material in the plane perpendicular to the nematic director and therefore represents a 
"discotic" phase. 



We first define the rescaled free energy Eq. ( 50 1 as 



9(S) 



127r^/?/(a0 
A3 



dt 
~t 



(1 + 0^ 



det[/ + t(/ + 7ct' 



+ c Tr cr' 



/-I 



Within the ansatz Eq. ( 56 1 it can be evaluated explicitly and is given by 



g{S) = -2 + 2 



v/3 + 7(l-^) 



' v/3 + 7(l-^) ^ 
^3 + 7(1 + 25) 



ln(i(3 + 7(l + 2^)))+^ 



6c 



(57) 



(58) 
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FIG. 2: The rescaled free energy difference (shifted for clarity by g{S — 0)) g{S) (Eq. (58l) for uniaxial ordering for the 
coupling parameter 7 — 0.5 for (from top to bottom) c — 0.000506 (isotropic), c = 0.000504 (coexistence between isotropic and 
nematic phases) c = 0.000502 (nematic phase). 



The first three terms represent the zero Matsubara frequency contribution to the electromagnetic field fluctuations. 
In case we would calculate the full Lifshitz expression in the case of magnetically inactive media, these three terms 
would need to be summed over all the Matsubara frequencies as in Eq. ^. 



We denote the first three terms in Eq. (58) as gvdw{S), i.e. 

9vdw{S) = g{S,c = 0). 

As one of the main aims of this paper is to investigate the effect of the full n-body van der Waals interactions it is 
interesting to compare this term with g).^ (5), which corresponds to the pairwise approximation, and can be derived 



from Eq. (58) or the expression Eq. ( |46[ ) truncated at the second order in a, which gives 

^?f.V(5) = i-^(5 + 4^')+0(7')- (59) 

The negative coefficient of S"^ clearly demonstrates that the van der Waals interactions favor ordering at the pairwise 
level. However, one notes that at the pairwise level there is obviously no difference between the energy of nematic 
ordering (ordering along an axis) and discotic ordering (concentration in a plane) because the interaction is even in 
S. In the study of Edwards and Gupta [T] it was found that within the ordered phase it is the nematic phase which 
has the lower free energy. The difference is due to the constraint part of g{S), corresponding to the last two terms of 



Eq. (58), denoted by gcoiS), and given by 



9co{S) = 9iS) - g^MS) = ^^^_^2S){1^- S) ' ^^^^ 

This term favors nematic ordering over discotic ordering as can be seen by its Taylor expansion to fourth order in S 

gco{S) = 9c + 18c5^ - 18c53 + 54cS'^ + ©(S^). (61) 

Here we can see how the symmetry of the van der Waals free energy between the nematic and discotic phases is 
broken by looking at the third order term in the interaction, corresponding to three body interactions. The total 
contribution to gvdw{S) to order 7^^ is then given by 

2 3 

5?iw(5) = i - + 4^2) + 2^(53 + 12052 ^ 10^3) ^ 0(^4)^ (02) 

Therefore, for small S, where the term dominates, the three body interactions tend to increase the van der Waals 
interaction energy and thus they reduce the tendency toward ordering. Furthermore, we also see that the three body 
interaction favors discotic rather than nematic ordering. 
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FIG. 3: The dependence of the c* and S* at the nematic transition on the interaction parameter 7 — ^ '^^'j^^^ j ■ 




FIG. 4: The dependence of the free energy difference g{S) ~ g{S = 0) on the order parameter 5* for various values of the volume 
fraction. Top to bottom v = 0.63,0.655,0.68. We have taken = 4, ex = 3 and ~ 80 and C = 2 x 10~®. The value for C 
corresponds to — 2.5 and A — 10. 



In Fig. tin we show the comparison of gvdw{S) with its approximation at third order in 7, gl^]\y{S), for the value 
7 = 0.5. We see that the approximation is very good at this value of 7, however the approximation tends to over 
estimate the zero Matsubara frequency van der Waals energy at all values of S. 

As far as the phase diagram based on the assumption of nematic ordering is concerned, the rescaled free energy 
g{S) for 7 = 0.5 is shown in Fig. ([2|. We see that when the value of c is sufficiently lowered a first order nematic 
phase transition is induced. The transition occurs at c* ~ 0.000504 and the order parameter S jumps discontinuously 
from to 5** ss 0.24. The complete dependence of the c* and S* on 7 is given on Fig. ([s]). 

The way in which c can be varied while keeping 7 constant is by changing the bare persistence length Ip of the 
polymer. The isotropic nematic transition occurs at a critical value I* = Sv/Giir'^X^c* . As I* is further increased, the 
order parameter S increases and eventually tends to 1. 

Inspite of the insights gained by this partitioning of the volume fraction dependence of 7 and c, it is necessary to 
investigate also the complete dependence of the free energy on the polymer density or equivalently its volume fraction. 
This dependence is hidden in 7 = 7(1^), Eq. (51), and c = c{i>), Eq. (53). 

Fig. ([4]) shows the order parameter dependence of the total free energy Eq. (58) for three different values of the 
volume fraction v contained in 7 = 7(1^) and c — c(j^). The dielectric constants were taken as ey = 4, = 3 and 
= 80, which corresponds to a polymer with a hydrophobic core in an aqueous solution. The persistence length 
was taken as 2.5 in the units of the monomer size and the radius of the polymer was taken as ten times the atomic 
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size cutoff, A = 10. 

For the chosen set of parameters the nematic transition occurs at v = 0.655, where the minimum of the free energy is 
displaced from its disordered value S" = to an ordered nematic state with 5^0. The critical concentration obviously 
depends on the exact values of the persistence length, the radius of the polymer and its dielectric anisotropy. 



VIII. DISCUSSION 



We explored the consequences of the orientationally dependent anisotropic vdW interactions on the ordering transi- 
tion of semi flexible polymers. A theory exact to all orders in the non-pairwise additive vdW interactions was derived 
and solved on the mean-field constant density level. We show that indeed there exists an orientational ordering 
transition engendering a nematic polymer phase. 

Through we limited ourselves to the zero Matsubara frequency term, the approach can be straightforwardly gen- 



eralized to include the full anisotropic Lifshitz expression [16]. In this case the free energy Eq. (581 needs to be 
rewritten as a full Matsubara frequency sum of the form 

9iS,i^)^X^ -2 + 2 :7^^=5 COS ( , ) + ^'M i^*-^ ^''^ + 

n=0 \ 

where now 



VV3 + 7n(l + 25) ; V3' "J J (l + 25)(l-5)' 

(63) 

( ^ ^ _ ^ (g||«n) -£-L««)) _ StTV _ 3 V 

In - - ^ ^ (^11 ^^^^^ _ ^^^^^^^^^ and c - ^ ^^^^^^^ - ^^^^^^ . [m) 

Depending on the dielectric anisotropy ej[(i^„) — £±{i£.n) at every Matsubara frequency, each term in the sum makes 
a contribution to the total ordering free energy. This means that the vdW part of the free energy in general increases 
and therefore displaces the ordering transition towards smaller densities. One should note that the overall density 



dependence of Eq. ( 63 1 becomes quite complicated in this case as partly (the first three terms) it depends on the 
frequency dependent dielectric response, and partly it doesn't (the last term). We also notice that the temperature 
dependence of the free energy upon including the non-zero Matsubara frequencies becomes much more complex. 

The total free energy thus contains two separate contributions: the Matsubara sum corresponding to the anisotropic 
vdW interactions and the last term stemming from the " constraint" free energy, that connects the dielectric anisotropy 
of the polymer chain with its configuration in space. It is the last term that can not be evaluated precisely because 
it contains a spatial cutoff A due to the continuum nature of the electrostatic field description. While this is not a 
major impediment regarding the scaling properties and the existence of the transition, it can introduce afterthoughts 
regarding the exact numerical values for the transition density. Based on numerical calculations we conclude that the 
polymer does not have to be very stiff in order to show the vdW driven orientational ordering transition. 

The mean-field approach outlined here for the case of orientationally dependent van der Waals interactions can 
of course be easily extended and applied to other cases. Mean-field electrostatic interactions would be a possible 
candidate. 
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